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ABSTRACT 

The gravitational waves (GWs) emitted by inspiraling binary black holes, expected 
to be detected by the Laser Interferometer Space Antenna (LISA), could be used to 
determine the luminosity distance to these sources with the unprecedented precision 
of <!%■ We study cosmological parameter constraints from such standard sirens, 
in the presence of gravitational lensing by large-scale structure. Lensing introduces 
magnification with a probability distribution function (PDF) whose shape is highly 
skewed and depends on cosmological parameters. We use Monte-Carlo simulations to 
generate mock samples of standard sirens, including a small intrinsic scatter, as well 
as the additional, larger scatter from lensing, in their inferred distances. We derive 
constraints on cosmological parameters, by simultaneously fitting the mean and the 
distribution of the residuals on the distance vs redshift (d^ — z) Hubble diagram. We 
find that for standard sirens at redshift z « 1, the sensitivity to a single cosmological 
parameter, such as the matter density Qm, or the dark energy equation of state w, is 
~ 50% — 80% tighter when the skewed lensing PDF is used, compared to the sensitivity 
derived from a Gaussian PDF with the same variance. When these two parameters 
are constrained simultaneously, the skewness yields a further enhanced improvement 
(by ~ 120%), owing to the correlation between the parameters. The sensitivity to 
the amplitude of the matter power spectrum, erg from the cosmological dependence 
of the PDF alone, however, is ^ 20% worse than that from the Gaussian PDF. The 
improvements for and w arise purely from the non-Gaussian shape of the lensing 
PDF; the dependence of the PDF on these parameters does not improve constraints 
relative to those available from the mean — z relation. At higher redshifts, the 
PDF resembles a Gaussian more closely, and the effects of the skewness become less 
prominent. These results highlight the importance of obtaining an accurate and reliable 
PDF of the lensing convergence, in order to realize the full potential of standard sirens 
as cosmological probes. 

Key words: cosmological parameters - cosmology: theory - gravitational lensing - 
gravitational waves 



1 INTRODUCTION 

The proposed space-based GW detector LISA, sensitive to 
frequencies between ~ 10~^ to ~ 0.1 Hz, will be able to 
detect massive black hole binary (MBHB) mergers out to 
redshifts «^5. In addition to providing information on black 
hole physics and general relativity, observations of GWs by 
LISA could be used as a probe of cosmology. As pointed out 
in a pioneering paper by Schutz (1986), GW observations of 
a binary system could yield an accurate estimate of the lu- 
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minosity distance to the source, independent of assumptions 
about the masses and orbital parameters of the binary mem- 
bers. Recent analyses in the context of LISA show that for 
many binaries, the luminosity distance could be measured 
to percent-level precision (see Arun et al. 2009b for a review 
and references). 

Although a typical LISA source will be relatively poorly 
localized on the sky (to ~ O.ldeg^), when the spatial infor- 
mation along the line of sight is taken into account, the num- 
ber of galaxies in the three-dimensional error volume will 
be reduced significantly (Holz & Hughes 2005; Kocsis et al. 
2006). Combining this with possible tell-tale time- variable 
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signatures (see Haiman et al. 2009 for a review of several 
possibilities), it may become feasible to identify an electro- 
magnetic (EM) counterpart. By measuring the redshift of 
the counterpart, a gravitational version of the Hubble di- 
agram could be constructed. This Hubble diagram, with a 
relatively small intrinsic scatter, and spanning a large range 
in redshift, can possibly impose tight constraints on cosmo- 
logical parameters (see, e.g., Arun et al. 2009a for a recent 
review and references). 

Unfortunately, however, gravitational lensing signifi- 
cantly changes this picture. Standard sirens, just as type la 
supernovae (SNe la), are (de)magnified by inhomogeneities 
in the matter distribution in the foreground, which intro- 
duces an uncertainty in the measured distance-redshift re- 
lation by up to ~ 10% for high-redshift (z^2) events (e.g. 
Holz & Hughes 2005; Kocsis et al. 2006). In the case of SNe 
la, proposed missions, such as the Supernova/Acceleration 
Probe (SNAP), are expected to find a few thousand use- 
ful sources. The random lensing magnification errors then 
average out, and even if they are unaccounted for, they 
have a relatively modest (^1%) impact on cosmological 
parameter-estimation, which can be ignored (Holz & Linder 
2005; Sarkar et al. 2008). The same strategy is unlikely for 
standard sirens: although the expected LISA MBHB event 
rate is highly uncertain, most models predict that it is sig- 
nificantly below the SNe la rate, with perhaps tens of de- 
tections per year (e.g. Menou et al. 2001; Sesana et al. 2007; 
Lippai et al. 2009; Arun et al. 2009b). Furthermore, the EM 
counterpart may be identifiable for only a fraction of these 
event so 

Motivated by the tremendous potential, in the absence 
of lensing, of standard sirens for cosmology, there have been 
proposals to correct for the effects of lensing of individ- 
ual sources on a case- by-case basis. These proposals in- 
clude measuring or constraining the magnification using ei- 
ther photometric and spectroscopic properties of foreground 
galaxies (e.g. Jonsson et al. 2007), or the combination of 
arcminute-scale shear and flexion maps (Shapiro et al. 2009; 
the earlier work of Dalai et al. 2003, which only consid- 
ered the shear, concluded that only modest, ~ 20%, cor- 
rections were feasible). Both of these methods could reduce 
the lensing-induced distance errors, in idealized cases, by a 
factor of up to ~two. 

Alternatively, several authors have investigated the pos- 
sibility of using lensing as a signal, rather than as noise, in 
probing cosmology (Dodelson & Vallinotto 2006, hereafter 
DV06; Linder 2008; see also Wang et al. 2009). In particular, 
the variance of the lensing magnification probability distri- 
bution function (PDF) depends on the amplitude of density 
fiuctuations and on the growth function, and therefore could 
be used to constrain cosmological parameters such as ug and 
Qrn- More specifically, DV06 showed that erg could be con- 
strained to an accuracy of ^ 5% by observations of 2000 

^ The Big Bang Observer (BBO), a concept for a space mission 
to succeed LISA, could detect a more than sufficient number of 
compact stellar binaries; this would even allow a useful measure- 
ment of the spatial power spectrum of the lensing convergence, 
which can provide additional cosmological constraints (Cutler & 
Holz 2009; see also Cooray et al. 2006 for the same idea with Type 
la SNe). In this paper, we will not consider the possibility of such 
a large number (^10'') of detectable events. 



SNe la. The overall shape of the lensing convergence dis- 
tribution contains additional cosmological information, be- 
yond the variance (Wang et al. 2009). For example, Linder 
(2008) emphasized that the theoretically possible minimum 
(de) magnification, which occurs along an empty beam, de- 
pends on cosmology. Several authors have indeed proposed 
fitting formulae for the lensing PDF, derived from numerical 
simulations, which are self-similar, and depend on cosmology 
only through the variance and the minimum magnification 
(see below). 

In this paper, we study the cosmological parameter con- 
straints from standard sirens, in the presence of lensing mag- 
nification. In general, lensing causes significant degradation 
of the constraints, which includes increased uncertainties, 
as well as a possible bias, in the parameters inferred from 
a given finite source sample. If the lensing PDF was Gaus- 
sian, and did not depend on the cosmological parameters, 
this degradation could be estimated simply by the increase 
in the variance of the inferred distance error to individual 
sources. However, the lensing PDF is highly skewed, and it 
does depend on the cosmological parameters. Our focus here 
is to quantify the extent to which these two features affect 
(and hopefully, mitigate) the degradation of the constraints 
expected from LISA. 

In the context of SNe, these questions have already 
been addressed in detail by several authors, including the 
impact of lensing on inferred dark energy parameters (e.g. 
Holz & Linder 2005; Sarkar et al. 2008) and on the normal- 
ization of the matter power-spectrum, erg (e.g. DV06). Here 
we consider, instead, a relatively small sample (~tens) of 
sources, with otherwise very small {<A%) distance errors, 
as expected from LISA standard sirens. The conclusions are 
not necessarily the same for such a standard siren sample as 
for the SNe la. This is because the probability distribution 
of the inferred distances, which is the convolution of intrinsic 
and lensing probability distributions, is much more skewed 
in the case of standard sirens, due to the dominance of the 
highly skewed lensing distribution. Additionally, depending 
on the statistic being used, having fewer events can increase 
the impact of non-Gaussianity on the inferred parameters. 

A few works have also studied the cosmological utility 
of standard sirens, and have included the effect of lensing, 
noting the significant degradation they cause in the con- 
straints (e.g. Holz & Hughes 2005; Dalai et al. 2006; Linder 
2008). Our present study adds to these existing papers in 
the following ways: (i) in our analysis, we include the de- 
pendence of the convergence PDF on cosmological param- 
eters, thus treating lensing as a potential source of signal, 
rather than as pure noise; (ii) we include erg among our pa- 
rameters, since the lensing PDF is particularly sensitive to 
this parameter; (iii) we perform multiple Monte Carlo real- 
izations of mock standard-siren samples, incorporating the 
non-Gaussian shape of the PDF, to obtain accurate esti- 
mates of the confidence intervals on the inferred parameters; 
and (iv) we use more recent fitting formulae for the lensing 
PDF, which improve the fit to numerical simulations. 

The remainder of this paper is organized as follows. 
In § [2j we briefiy summarize the basic background mate- 
rial required for our study, including information on both 
standard siren distance measurements and on gravitational 
lensing. § [3] outlines the details of our Monte-Carlo simu- 
lation procedure. In § U we present and discuss our main 
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results, and contrast these with the analogous results in 
the case of SNe la. Finally, we summarize our main con- 
clusions in § O Throughout this paper, we adopt stan- 
dard ACDM as our fiducial cosmological model, with pa- 
rameter values consistent with the WMAP fifth year re- 
sults (Komatsu et al. 2009), i.e., {Q,m,^A,^b,h,<js,,ns} = 
{0.279,0.721,0.0462,0.701,0.817,0.96}. These values are in 
general agreement with most recent, seventh year results, as 
weU (Komatsu et al. 2010). 



2 STANDARD SIREN DISTANCE 

MEASUREMENT AND GRAVITATIONAL 
LENSING 

We follow the convention in the astronomical literature, and 
express the luminosity distance dL{z), inferred from an ob- 
servation of a standard siren at redshift z, by the distance 
modulus 



value, we adopt the same strategy as in Kocsis et al. (2006), 



m,{z) = mo{z) + 5mint(z) + (Smicn(z). 



(1) 



Here mo{z) is the true distance modulus in a homogeneous 
universe, depending only on the usual luminosity distance 



mo{z) = 51ogio 



dLjz) 
10 pc 



(2) 



and in a fiat universe as assumed throughout this study, 
dz' 



dL{z) = {1 + z)c 



Hiz')' 



(3) 



where c is the speed of light, and H{z) is the Hubble pa- 
rameter. 

The other two terms in equation Q vary from source to 
source, and account for the intrinsic dispersion due to LISA's 
instrumental and background confusion noise (5mint), and 
for the scatter caused by lensing magnification (Sniicn)- In 
equation an arbitrary calibration constant has been 
omitted. We furthermore ignore contributions to the scat- 
ter from peculiar velocities, which are only important at 
low redshifts (z_^0.3) and also neglect any other source of 
systematic errors, which, of course, need to be included in 
the actual data analysis. 

For a standard siren, the intrinsic error Smint(z) is usu- 
ally very small, but the distribution of 5mint{z) depends on 
properties of the detector and of the MBHB in a non-trivial 
manner. In the most general case, the width of the intrinsic 
error distribution is a function of 17 parameters (Vecchio 
2004; Holz & Hughes 2005; Kocsis et al. 2006; Trias and 
Sintes 2008). Owing to the complexity of the problem and 
to computational limitations, the distribution of intrinsic lu- 
minosity distance errors has been estimated for only a few 
discrete choices of redshifts and BH masses. To scale the in- 
trinsic error to an arbitrary redshift and mass from a base 



We note that in some applications, using fluxes, rather than dis- 
tance moduli, is preferable, since lensing acts linearly on the flux, 
and therefore it does not shift the mean of a flux PDF (whereas 
lensing can shift the mean of the magnitude PDF). However, as 
explained below, in our analysis, we fit the entire PDF, and the 
two prescriptions would lead to the same result. 



SdL{M,z) 

dL{M, z) 



{S/N){M,z) 



{S/N){Mo,zo) 



SdL{Mo,zo) 
dL{Mo,zo) 



(4) 



where S/N is the expected signal-to-noise ratio (SNR) for 
the detection of GWs. Equation (|4|) is motivated by the fact 
that the signal amplitude is inversely proportional to the lu- 
minosity distance. In computing the SNR, we have assumed 
one year observation time prior to coalescing, and used the 
same noise spectrum as in Lang & Hughes (2006), which 
includes instrumental nois^fl and confusion noise (from both 
extragalactic and galactic sources). Limited by the low- 
frequency noise wall, the visible time of very massive MB- 
HBs can be less than one year, rendering equation Q in- 
accurate at large BH mass M. Nevertheless, this, as well as 
other details of the intrinsic error distribution, make very 
little difference to the results of the present study, as the in- 
trinsic error is smaller than the lensing error. The base value 
of the luminosity distance error is chosen to be 2.88 x 10"'^ 
for an equal- mass binary (Mi = M2 = 3 x 10^ M0) at z = 1, 
computed by Klein et al. (2009) using the full second post- 
Newtonian (2PN) gravitational waveform, and also includ- 
ing spin-orbit precession. The variance of Srriint, cfnti is re- 
lated to the fractional distance error SdL/dL through equa- 
tion (O, 



(Tint = 551ogiodL « 2.17 



Sdj_ 

dL 



(5) 



The lensing error 5micn(z) is directly related to the lens- 
ing magnification /i. 



(5mi, 



~2.51ogio^ = 51ogio(l 



(6) 



where k is the convergence. The PDF of P{k,), is known to 
be skewed, since there exists a minimum convergence, Kmin, 
corresponding to an empty path between the source and the 
observer. Since in this paper, we are primarily interested 
in quantifying the effects of the non-Gaussian tails of the 
PDF, we have to rely on estimates of the PDF in cosmolog- 
ical simulations. Using the results of N-body simulations in 
White (2005), Wang et al. (2009) compared three different 
models for P{k), proposed by Taruya et al. (2002, i.e. the 
log-normal distribution), Wang et al. (2002, hereafter W02) 
and Das & Ostriker (2006, hereafter DO06). They concluded 
that the fitting formulae proposed by DO06, Poooe best fit 
the simulation data, while the formulae proposed by W02, 
Avo2, over-predict the tail of the convergence distribution. 
In the following, we therefore adopt Pdoo6 as the "true" 
distribution, 

Pt,uo(K) = Pdoo6(k) = N (7) 
[\n{l + ^/\n^in) + 'S^/2\)f 



X exp 



2E2 



X [l+4/(l + K/|Kminl)]} 

X ■ . 

min I ^~ ^ 

Here, the three parameters A'^, E and A are fixed by three 



^ Generated by the on-line sensitivity curve generator, 
|http://www.srl. caltech.edu/~shane/sensitivity I 
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Figure 1. PDFs of lensing convergence for sources at 2 = 1 (up- 
per panel) and z = 4 (lower panel), from fitting formulae proposed 
by Das & Ostriker (2006) and by Wang et al. (2002), compared to 
a log-normal and a Gaussian distribution, as labeled. Note that 
i-Vv02 is sharply truncated in order to satisfy the constraint equa- 
tion lITOt : the cutoff oi z = 1(4) is at k = 0.25 (1.46, outside the 
K range of the figure). The low— redshift PDFs are more skewed 
than the high redshift ones, and Pvv02 predicts a significantly 
larger high-re tail than PdO06- 



constraint equations, 



K.Ptiuc{K.)dli — 0, 



^-Ptruc(«:)dK = 



(8) 

(9) 
(10) 



where is the variance of k. Note, in particular, that Pdoo6 
depends on cosmology only through Kmin and a^. These, in 
turn, are determined from the equations 



3„ fHoV ,^ 

Kmin = — -Sim I — I J ax(l 



and 

2 



9n fQmHi^' 



dx 



+ 



xixs - x) 



Xs 



(12) 



where x is the comoving distance to redshift z, Xs is the 
comoving distance to the source, and A^(fc, z) is the dimen- 
sionless non-linear matter power spectrum, evaluated in this 
study using the algorithm of Smith et al. (2003). Similar to 
equation ([5|, the variance of Srriicn could then be computed 
from (j^. In the case of a Gaussian distribution, using equa- 
tion ((6)1, we have 
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(13) 



LU 

0) 




Figure 2. The standard deviation of the intrinsic and lensing 
error distributions. The latter assumes an equal-mass binary 
(Ml = M2 = 3 X IO^Mq). In contrast to SNe la, the intrinsic 
error for a typical standard siren at redshift 2<^4 is a factor of 
^two smaller than its lensing error (and nearly an order of mag- 
nitude smaller at z ~ 1). 



It is important to note here that, in comparing fits to 
simulation data, Wang et al. (2009) treated Kmin as a free 
parameter. When Kmin is fixed at its theoretical value from 
equation (|lip . Pdoo6 was, in fact, found to under- predict 
the high-K tail of the PDF. The true PDF then lies in- 
between Pdoo6 and Pwo2. Besides the high-K tail, other 
issues with Pwo2 include discontinuities and violating flux 
conservation at (Zentner & Bhattacharya 2009). We 

therefore adopt the fitting formulae from DOG, despite their 
imperfect fit to the simulation data. To illustrate the dif- 
ferences between the various PDFs, in Figure [1] we show 
Pdoo2(k), Pwo2(k), a log- normal PDF, and a Gaussian PDF 
for sources at z = 1 (upper panel) and 2 = 4 (lower panel). 
The cutoff of Pwo2(k) at k = 0.25 is imposed in order to 
satisfy the constraint equation (|10p : Pwo2(k) of 2; = 4 is 
also truncated but the truncation is outside the k range 
shown in the figure. As Figure [1] shows, Pdo6 and Pwo2 
are both skewed, and the skewness is more pronounced at 
lower redshift. In Figure [2j we show the r.m.s. magnitude 
of the intrinsic and lensing dispersions as a function of red- 
shift. For the latter, we assumed an equal-mass binary with 
Ml — M2 = 3 X IO^Mq, which is close to the best case for 
LISA. At 2:<;4, the intrinsic dispersion is always a factor of 
^two smaller than the lensing dispersion. 

In reality, there will most likely be a wide distribution 
of masses and mass ratios among the BH binaries detected 
by LISA. These distributions are highly uncertain, and for 
simplicity, we avoid modeling it in this paper. When both 
masses are in the range of ~ 10^ — 10® , the variation of 
the intrinsic error with the mass ratio is relatively modest 
(within a factor of ~ 2; Klein et al. 2009), and the fiducial ac- 
curacies we chose are typical of these binaries. For binaries 
with component masses significantly outside this optimal 
range, the intrinsic error will be larger than assumed here, 
and for some events, it will therefore likely become compara- 
ble (or larger) than the lensing error. Nevertheless, for most 
events, we expect the lensing error to remain dominant, at 
least at redshifts as low as 2; ~ 1, where the differences shown 
in Figure [2] are nearly an order of magnitude. 
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Figure 3. Probability distribution of the total magnitude error of 
standard candles (solid [red] curve) and of standard sirens (dashed 
[green] curve) at z = 1 . The intrinsic errors of both standard can- 
dles and of standard sirens are assumed to have Gaussian distri- 
butions (in magnitude), with a standard deviations of 0.1 mag, 
and 2.88 X 10~^ mag, respectively; this Gaussian is convolved 
with the DO06 lensing PDF (shown separately in Figure [TJ. 

In Figure [3l we show the overall probability distribu- 
tion of the magnitude error, given by the convolution of the 
distributions of Srriint and Smien- The solid (red) curve is 
for standard candles, assuming an intrinsic dispersion that 
is Gaussian (in magnitude) with a standard deviation of 0.1 
mag (a canonical value adopted by DV06 and in many other 
papers; e.g. Holz & Linder 2005; Zentner et al. 2009). The 
dashed (green) curve is for standard sirens, again assum- 
ing a Gaussian intrinsic error distribution with a standard 
deviation of 2.88 x 10~^ mag. In both cases, the source is 
assumed to be at redshift z = 1. As the figure shows, the 
overall error distribution of standard sirens remains much 
more skewed than that of standard candles. 

We note that the intrinsic distribution of Srriint may well 
be skewed, as well. However, as long as Smint is several times 
smaller than Smian, the non-Gaussianity in Sm-mt should not 
have a significant impact on our conclusions. 



3 LIKELIHOOD ANALYSIS AND 
MONTE-CARLO SIMULATION 

We next use Monte-Carlo (MC) simulations to derive confi- 
dence intervals on the parameters inferred from mock sam- 
ples of standard sirens. We follow a procedure similar to 
DV06, consisting of the following steps: 

(1) First, we generate random redshifts for Ncvant stan- 
dard sirens, distributed uniformly in a bin of width Az = 1, 
centered at redshift Zccnt- The number of events, A'ovcnt, 
in each realization is varied, ranging from 40 to 2000. The 
lower end of this range corresponds approximately to a (pes- 
simistic) expectation for the total number of LISA events, 
while the upper end corresponds to the number of SNe ex- 
pected in future SN surveys such as SNAP/JDEM (allowing 
us to contrast standards sirens to SNe). The redshift 2:ccnt 
is set at 1, 2, 3, and 4, i.e., the events are between redshift 
0.5 to 4.5. At z<^0.5, very few standard sirens are expected, 
and the errors from peculiar velocities become comparable 
to those from lensing (e.g. Kocsis et al. 2006) and would 



have to be included in the analysis. At z^4.5, the lensing 
PDF approaches a Gaussian shape, and the effects of non- 
Gaussianity become negligible, as we will see in § [J] 

(2) For each event i, we compute the true distance mod- 
ulus using equation We then draw two random values of 
K, one from the distribution of Ptruc(fi:) and the other from 
f'Gaaa('v) whosB Variance is set to be the same as predicted 
in equation (jlOp . and convert these k values to Smi^n- We 
also draw a random Srriint from a Gaussian distribution with 
zero mean and variance (Tint- 

(3) We add Srriint and Smicn to mo to form two sets of 
mock data, {zi, mi), corresponding to the two choices of the 
convergence PDF. 

(4) We then perform two separate likelihood analyses 
on the mock data, assuming the convergence PDF, either 
true or Gaussian, is known. The distribution of Srriint is then 
convolved with that of Sniicn to obtain the PDF of the total 
error, Ptot. In this analysis, in difference from DV06, we as- 
sume that (Tint is known and is fixed at the "correct" value 
(i.e. the same value used in step (2) above to generate the 
mock sample) . As mentioned above, in our case, the intrinsic 
dispersion is sub-dominant compared to the lensing disper- 
sion, and we therefore do not expect marginalizing over this 
intrinsic dispersion will modify our results. 

(5) In each analysis, we vary either one or two cosmo- 
logical parameters, and find the parameter value(s), p™^", 
that maximizes the likelihood C, or equivalently minimizes 
the quantity H = — 21n£. For a normalized Ptot, 

H = -21nZ: = -2^1nPtot(m, -mo{z,)), (14) 

i 

where terms independent of cosmology have been omitted 
(Marshall et al. 1983). In the limit of Gaussian probability 
distribution Ptot, —2lnJZ becomes similar to the usual 
statistic. 

6) We repeat the steps (l)-(5) 1000 times to produce two 
distributions of p™'"', corresponding to the two distributions 
of K. 

7) Finally, we evaluate the effects of the convergence 
PDF on the cosmological constraints by comparing these 
two distributions of p™^" . 

The above approach allows us to quantify the impact 
of the non-Gaussian shape of the convergence PDF on the 
distribution of the parameter estimates. It is useful to con- 
trast here the above approach to some others that are of- 
ten used in the literature to forecast parameter errors. The 
Fisher-matrix formalism (Tegmark et al. 1997) is commonly 
adopted when the number of parameters to be estimated 
simultaneously is too large. This method assumes that the 
likelihood function is Gaussian, both as a function of the 
parameters to be estimated and as a function of the observ- 
ables. In the present case, both of these assumptions are 
false (our results will be explicitly compared to those from a 
Fisher matrix below). Another common practice is to com- 
pute the likelihood, as a function of the parameters, around 
a single realization of the mock data - usually, either the 
most likely realization (see, e.g. Figure 6 in Holz & Hughes 
2005), or in a "typical" realization (see, e.g., DV06). This 
avoids the assumption that the likelihood is Gaussian in the 
parameters, but does not fully capture the non-Gaussian 
shape of the distribution as a function of the observables. 

In the context of SNe, Holz & Linder (2005) have used 
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full Monte Carlo simulations to quantify the impact of lens- 
ing on the distribution of the inferred best-fit parameters, 
similar to our procedure above. Because of the large intrinsic 
dispersion of the SNe, lensing is less important to begin with, 
and the effect of the non-Gaussianity is also relatively mod- 
est (see their Figures 8 and 9). As far as we are aware, the 
same approach of using full Monte Carlo simulations have 
not been used to study the impact of lensing for standard 
sirens, for which the intrinsic dispersion is smaller, and for 
which the effects of lensing, as well as of non-Gaussianities, 
are more important. 



4 RESULTS AND DISCUSSION 

4.1 Sensitivity to a single cosmological parameter 

We begin by investigating the sensitivity of standard sirens 
to a single parameter. We focus on the degradation caused 
by lensing, and how this degradation is affected (hopefully, 
mitigated) by the non-Gaussian shape of the lensing PDF 
and by its dependence on cosmological parameters. 



4-1.1 Matter density fim 

In the first example, we study the constraints on Clm while 
other parameters are held fixed at their fiducial values. The 
number of events in each realization is 40, and Zcont is set 
either at 1, 2, 3, or 4. The distributions of best-fit Qm val- 
ues are shown in Figured) where the solid (red) curves and 
dashed (blue) curves correspond to the results when Ptrue(K) 
and PGaaa(K) are used in the likelihood analysis, respectively. 
The shaded (green) region also shows the distribution for 
PtTue{ii), but with the PDF fixed using the fiducial values 
of the cosmological parameters. The true (fiducial) Qrn is 
indicated by the vertical dotted lines. 

The parameters of these distributions are summarized 
quantitatively in Table ([1]), which lists the peak, 68% and 
95% errors of best-fit flm and the ratios of the errors re- 
sulting from the two different assumptions about the con- 
vergence PDF (Gaussian vs. non-Gaussian). The second to 
last column shows the degradation due to the lensing, indi- 
cated by the ratio of 68% errors with and without lensing. 
For a comparison of constraining power of individual stan- 
dard sirens and that of SNe la, the last column of Table 
IT} shows how many standard sirens need to be observed in 
each redshift bin in order to obtain the same constraints on 
flm as two thousand SNe la at redshifts between 0.5 and 
1.5. Within the parentheses are the corresponding numbers 
of events but in the hypothetical case that lensing is either 
absent or can be perfectly corrected for. In particular, an 
event in the first redshift bin on average has cosmological 
constraining power similar to 2000/114 ~ 18 SNe la. 

Several interesting conclusions can be drawn from Fig- 
ure|4]and Table[T] First, consistent with previous works, our 
results show that lensing severely degrades the constraints 
from standard sirens, and that this degradation is worse than 
for SNe la, due to the small intrinsic dispersion of standard 
sirens. At 2 ~ 1, lensing increases the 68% error by a factor 
of 4.71 for standard sirens, while only by 25% for SNe la. 
We also note that the best-fit values of i},n are essentially 
unbiased; this is simply because we have assumed that the 



shape of the lensing PDF is known ab-initio (we will relax 
this assumption in § 14.41 below) . 

Second, the constraints on i},n become worse as Zcont 
increases. This is because both the intrinsic and the lensing 
dispersions rapidly increase with redshift, as shown in Figure 
([21 . When the number of events in a realization is fixed, the 
error on Qrn is roughly proportional to the total dispersion 
of the measured distance modulus, 



Sm 
2fl, 



dm 



dQrr. 

5m. 



oc 5m 



1 dd, 



dL dQ,r, 



(15) 



In the above, we have used the approximation H{z) « 
//o\/f2m(l -|- z)'^ (z^l), which allows us to move Q,m out of 
the integral of d_L, and cancel the redshift dependence. Tak- 
ing the values of crtruc in Table ([l]), and using the scaling that 
the error on a cosmological parameter is inversely propor- 
tional to the square root of the number of events, our results 
show that an event at z — 1 has roughly the same contribu- 
tion in constraining fi^ as (1.14 x 10"^/3.49 x 10"^)^ « 11 
events at z = 4. 

The third and most interesting conclusion is that the 
constraints in the true PDF case are tighter than those in 
the Gaussian PDF case. The difference is as large as ~ 50% 
in the lowest redshift bin, where the true PDF is the most 
skewed and is furthest from a Gaussian shape, as clearly seen 
in Figure [1] This illustrates the importance of the shape of 
the convergence PDF, which, if not appropriately treated, 
could result in a biased estimate of the constraints. For 
example, the popular Fisher matrix technique, as we have 
checked, gives constraints on Qm that are very close to our 
Gaussian case. For any attempts of correcting for lensing on 
a case-by-case basis, this means that not only the variance, 
but the full shape of the PDF needs to be considered in the 
analysis, in order to evaluate the gain from such corrections. 
In general, these results are good news for the prospects of 
using standard sirens in cosmology, as the degradation from 
lensing is somewhat less severe than estimated from the vari- 
ance alone. We also note, however, that the benefits of the 
non-Gaussianities largely disappear at redshifts beyond z^2, 
where the lensing PDF is closer to a Gaussian. 



4-1.2 Dark energy equation of state w 

We have also run MC simulations in which we varied the 
dark energy equation of state parameter w, while holding 
all other parameters fixed. The results are shown in Table 
([21 and Figure [21 The mock survey parameters are the same 
as in Table [1] The results are qualitatively similar to those 
in the case of Qrn- The benefits of the non-Gaussian shape 
are, however, somewhat larger, with the degradation due to 
lensing smaller by up to ~ 80% (in the lowest redshift bin) 
compared to the Gaussian case. We also note that the errors 
on w have a steeper redshift dependence than those of Qrn- 
An event ai z — 1 has a constraining power on w equivalent 
to (6.83xl0~^/8.08xl0"^)^ ^ 71 events at z = 4, compared 
to ~ 11 for Qrn- The transition from cosmic deceleration to 
acceleration happens at 2 ~ 1; events at this redshift are 
therefore especially sensitive probes of dark energy. 
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0.28 0.29 

Figure 4. Distributions of Q,m that maximize the likelihood. The solid [red] curve, and the dashed [blue] curves are the distributions 
when Ptruc(«;) and PcassC'^) are used, respectively. The shaded [green] region also shows the distribution for Ptruo(K), but with the PDF 
fixed using the fiducial values of the cosmological parameters. The vertical dotted lines in this and following similar figures indicate the 
fiducial value of the cosmological parameter being constrained. The number of events in each realization is 40. 



Table 1. The peak position (y) and the 68% and 95% confidence levels for the distribution of the best-fit f2m when the convergence PDF 
is assumed to be Gaussian (denoted by the subscript "Gass") and when it is assumed to follow the non— Gaussian DO06 distribution 
(denoted by the subscript "true"). The ratios of the 68% errors in these two cases is also shown. Each realization generates a mock 
sample of 40 standard sirens, distributed uniformly in redshift within a width Az = 1, centered at Zcent = 1, 2, 3, or 4. The second to 
last column shows the degradation caused by lensing by comparing the 68% errors with and without lensing effects. The last column 
lists the number of standard sirens that would have similar constraining power (evaluated using the 68% error) as 2000 SNe la uniformly 
distributed between redshift 0.5 and 1.5. Within the parentheses is the corresponding number of events in the absence of lensing or with 
the lensing effects perfectly corrected. 



•Scent 




Gaussian PDF 

'''68, Gass 0'95,Gass 






True PDF 

0"68,truc 


fgS.true 


'^68,Gass 


0-95,truc 


<^68,nolon5 




1 


0.279 


5.22 X 10-=* 


1.06 X 10" 


-2 


0.279 


3.49 X 10^3 


7.11 X 10-3 


1.50 


1.49 


4.71 


1.14 X 102(8) 


2 


0.279 


7.06 X 10-^ 


1.42 X 10" 


-2 


0.279 


6.24 X IQ-^ 


1.26 X 10-2 


1.13 


1.12 


4.57 


3.65 X 102(27) 


3 


0.279 


8.60 X 10-3 


1.73 X 10- 


-2 


0.279 


8.52 X 10-3 


1.74 X 10-2 


1.01 


0.99 


3.00 


6.80 X 102(118) 


4 


0.285 


1.14 X 10-2 


2.27 X 10" 


-2 


0.276 


1.14 X 10-2 


2.31 X 10-2 


1.00 


0.98 


1.83 


1.21 X 103(565) 



4-1.3 Power spectrum normalization ag 

We also studied as, which is different from Qm and w in that 
as has no effect on the mean distance modulus, and therefore 
could not be as tightly constrained by LISA alone. To obtain 
good constraints (for numerical convenience), four hundred 



events are generated in each realization, and the results are 
shown in Figure [6] and Table ((3| . Note that constraints on 
as, as opposed to those on Qm and w, have only a very 
mild dependence on redshift: the 68% errors are ~ 0.025 
at z = 1,2 and 3 for the case of the Gaussian distribution. 
This result could be understood by the following argument. 
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Figure 5. Distributions of best-fit w values when the true (solid [red] curve) and Gaussian (dashed [blue] curve) PDF are used in the 
likelihood analysis. 1000 mock samples are generated with 40 events in each realization, distributed randomly in redshift within a bin 
centered at Zccnt = 1, 2, 3 or 4, as labeled. As in Figure|4] the shaded [green] region shows results when the lensing PDF artificially uses 
the fixed values of the fiducial cosmological parameters. 

t 

Table 2. Peak, 68% and 98% errors of best— fit w values, for different assumptions about the convergence PDF, as in Table [T] The ratios 
of the errors for the different PDFs are also shown. Each realization includes 40 events. 
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^truc 


True PDF 
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""95, true 




CT68,truc 


0-95,Gaas 
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-1.000 


1.43 X 10-2 


2.94 X 10" 
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-1.000 


8.08 X 10-=^ 


1.66 X 10" 
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1.78 


1.77 
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-0.994 


3.01 X 10-2 


6.06 X 10" 


-2 


-1.006 


2.33 X 10-2 


4.73 X 10" 


2 


1.29 


1.28 


3 


-1.000 


4.76 X 10-2 


9.67 X 10- 


-2 


-1.000 


4.31 X 10-2 


8.85 X 10- 


2 


1.11 


1.09 
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-0.964 


7.32 X 10^2 


1.47 X 10" 


1 


-1.024 


6.83 X 10-2 


1.38 X 10- 


1 


1.07 


1.06 



First, we see from equation (|12|) that oc A'^{k, 2) oc of 
which the latter proportionality is only approximate because 
of nonlinear effects. This yields 

^-4=^-4= ^^'^--f''' = (2/iV_0^^' - 0.071, (16) 

independent of redshift. When estimated from equation (|16p 
instead, 5as ~ 0.0236, which is only slightly lower 

than the values from the numerical calculation; the small 
difference can be attributed to the intrinsic dispersion. 

Another important difference is that non-Gaussianities 



degrade, rather than mitigate the lensing errors on erg, com- 
pared to the Gaussian PDF case - by ~ 20% in the lowest 
redshift bin. This is in contrast to the situation for Q.m and 
w, for which non-Gaussianities tighten the constraints. The 
reason for this finding will be investigated in § 14.31 below. 



4.2 Information from the parameter-dependence 
of the PDF 

As explained above, we have included two effects about the 
lensing PDF in our analysis: (i) its non-Gaussian shape, and 
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Figure 6. Distributions of best-fit erg values in 1000 Monte-Carlo realizations when the true (solid [red] curve) and a Gaussian (dashed 
[blue] curve) PDF are used in the likelihood analysis. 400 events are generated in each realization, in a redshift bin centered at ^cent = 1, 
2, 3, or 4, as labeled. 



Table 3. Same as Table l[2]l, except that eg is being constrained and each realization produces 400 events. 
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0.814 


2.42 X 10-2 
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2 


0.814 


2.93 X 10-2 


5.76 X 10" 
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0.814 


2.48 X 10-2 


4.76 X 10" 
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0.814 


2.78 X 10-2 


5.32 X 10" 
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0.89 


0.89 
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0.814 


2.64 X 10-2 


5.21 X 10" 
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0.814 


2.83 X 10-2 


5.74 X 10~ 
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0.93 


0.91 


4 


0.826 


3.56 X 10-2 


7.01 X 10" 
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0.814 


3.78 X 10-2 


7.45 X 10" 


2 


0.94 


0.94 



(ii) its dependence on cosmological parameters. It is impor- 
tant to clarify which of these two effects was responsible 
for the improvement in the constraints on Q.m and w, rel- 
ative to the case when the lensing PDF is Gaussian and is 
considered pure noise (i.e. with no parameter-dependence). 
Assuming a Gaussian PDF and employing the Fisher ma- 
trix method, Zentner & Bhattacharya (2009), find that the 
dispersion of SNe la moduli can help break degeneracies be- 
tween dark energy parameters. However, the non-Gaussian, 
high-convergence tail of the lensing PDF has been shown 
explicitly to contain cosmological information, whose statis- 
tical accuracy in an all-sky map is competitive with other 
cosmological probes (Wang et al. 2009). 

We perform two tests in order to establish which of the 



two effects is more important. First, we artificially fix the 
shape of the PDF with the fiducial values of cosmological 
parameters throughout the analysis. The results of this aca- 
demic exercise are shown in Figures |4] and [5] by the shaded 
regions. This brings only a marginal change in the distri- 
butions of the best-fit values compared to the fiducial case, 
showing that the effects of the true lensing PDF come pri- 
marily through its non-Gaussianity, rather than through its 
cosmology-dependence. 

This conclusion is confirmed by a second test, in which 
we artificially fix the distance modulus, and re-derive con- 
straints on Q.m- We find that the constraints in this case 
become by about an order of magnitude worse. A closer in- 
vestigation shows that for Q.m, the cosmological dependence 
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of the PDF, in fact, slightly harms the constraints, while 
for w, it slightly improves the constraints. The fact that the 
cosmological dependence of the PDF harms the constraints 
on Qrn must also be caused by the non-Gaussian shape of 
the PDF. This is because for a Gaussian PDF, it could be 
argued (e.g. using the Fisher matrix) that additional infor- 
mation can only improve the sensitivity to a parameter (by 
increasing the absolute value of the corresponding element 
in the Fisher matrix) . We have verified by numerical calcula- 
tion that for a Gaussian PDF, the cosmological dependence 
indeed improves constraints on both f?™, and w, in agree- 
ment with this argument. 



4.3 Why does non-Gaussianity of the PDF affect 
the constraints? 

In this section, we discuss why the parameter constraints are 
affected by the shape of the convergence PDF. To illustrate 
this, we consider a quantity f , defined as, 

^ = = -2Ptot(5mtot)lnPt'ot(5m;„0, (17) 

where (H) is the expectation value of H for one event. In 
equation (|17|) . the factor Ptot(5mtot), computed with the 
true convergence PDF and with the fiducial values of the cos- 
mological parameters, is the probability that an event has a 
(true) total deviation of Snitot- The factor — 21nPt'ot('5"itot)) 
computed with the convergence PDF adopted for the like- 
lihood analysis (i.e. either a Gaussian or the true PDF) 
and with the varying cosmological parameters, represents 
the contribution to H from a single event with this Smtot- 
Note that the residual distance modulus inferred in a test 
cosmology, 5m[at, differs from the true residual 5mtot, since 
the luminosity distances in the two cosmologies are differ- 
ent. In the top left panel of Figure [71 we show ^ at z = 1 as 
a function of the true residual Smtot, where Qm is set at its 
fiducial value. 

The change in ^ when Qm is varied, A^, is shown in 
the left middle and left bottom panels, in which Qm is de- 
creased or increased by ues.Gass (= 5.22 x 10"^, see TableO, 
respectively. The main effect of changing Qm is to shift the 
measured residual Sm^^t away from its true value of Smtot, 
producing wave-like shapes in the curves — A^{Smtot)- 
For example, assuming an Qm larger than its fiducial value 
(bottom panel) increases the apparent residuals Sm't^t by re- 
ducing the distances (^^(z)) and therefore reducing mo. In 
other words, the source is inferred to be too close, and the 
magnification is inferred (statistically) to be smaller than 
the correct value. This corresponds to a reduced probability 
and increased ^ at large positive Smtot where the Ptot (Smtot) 
decreases with Smtot (i.e. if the source was already strongly 
demagnified; see the dashed curve in Figure [3|. Conversely, 
it results in an enhanced probability and a decreased ^ at 
negative Smtot (if the source was strongly magnified, it is 
inferred to need less magnification, which has a larger prob- 
ability). The increase of at mildly positive Smtot and the 
decrease at mildly negative Smtot are both stronger in the 
case of true PDF. This is because near Smtot ~ 0, the con- 
tribution to 5 is a much steeper function of Smtot for the true 
PDF, as seen clearly in the top panel. Indeed, the change of 
E, i.e. the integral of A^, is correspondingly also larger in 



the case of the true PDfQ, leading to an improved constraint 
on Qm- 

The fact that the sharp features of the non-Gaussian 
lensing PDF near its peak help tighten the sensitivity to 
may be surprising, since they are in contrast with the more 
widely discussed, detrimental, effects of long non-Gaussian 
tails. For example, the mean of A'^ numbers, generated from a 
non-Gaussian distribution, convergences less rapidly (than 
1/y/N) to the correct value, due to such shallow, non- 
Gaussian tails. However, when one fits the entire shape of 
the distribution (rather than utilizing just the mean), then, 
intuitively, the benefits of sharp features are not surprising. 
For example, one can envision a hypothetical, picket-fence- 
shaped distribution with the same r.m.s. as a smooth Gaus- 
sian. It is easy to see that if the picket-fence pattern shifts 
with a parameter, then this parameter can be arbitrarily 
tightly constrained, in the limit that the picket-fence shapes 
become arbitrarily sharp Dirac J-functions. 

The right panels of Figure [7] show the same quantities 
as in the left panel, except for era. Unlike Q,m, o"8 does not 
change the mean Sm't^f As a result, in the Gaussian case, 
the curves of A^ are strictly symmetric about Smtot = 0. 
The middle right and the bottom right panels show that the 
non-Gaussianity again enhances the sensitivity to as near 
the peak (for modest demagnifications) , whereas ^ in the tail 
of the true PDF becomes relatively insensitive to the change 
in (T8. Overall, this loss of sensitivity in the tails is the larger 
effect, and it renders AE smaller than the Gaussian cas^f]. 
This explains why the constraints on ag, becomes worse when 
the true PDF is used. 

4.4 How important is it to know the correct 
shape of the PDF? 

Our results above show that the shape of the lensing PDF 
has a non-negligible effect on the parameter constraints. 
However, the analysis above has assumed that the shape 
of the lensing PDF is precisely known. 

A related question of interest, therefore, is: how accu- 
rately do we need to know the lensing PDF, in order to avoid 
either a bias in the inferred parameters, or a mis-estimation 
of its confidence levels? 

In order to investigate this issue, we here create mock 
data-sets with the "true" lensing PDF from DO06, but we fit 
them using "wrong" PDFs. The two "wrong" PDFs we chose 
are a Gaussian (representing an extremely wrong case) , and 
a log-normal distribution (which represents a more reason- 
able estimate of the current uncertainty about the conver- 
gence PDF). The results for Q,m, w and as are shown in 
Tables HHg] 

Tables [4] shows that the maximum-likelihood estimate 
(MLE) of Q.m remains essentially unbiased: the MLE distri- 
bution peaks at values that are offset from the true value 

■* Specifically, for the cases shown in the left panels of Figure [T] 
we find AS= 0.0454 (true PDF, reduced), 0.0232 (Gaus- 
sian PDF, ilm reduced) ,0.0506 (true PDF, increased), 0.0214 
(Gaussian PDF, Vlrn increased). 

^ For the cases shown in the right panels of Figure [7] AS= 
0.0190 (true PDF, H™, reduced), 0.0289 (Gaussian PDF, re- 
duced), 0.0167 (true PDF, increased), 0.0259 (Gaussian PDF, 
Q,m increased). 
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Table 4. The peak of the maximum-likelihood estimator of and its 68% and 95% errors when the mock data, produced using the 
"true" PDF, is fitted by erroneously using either a Gaussian or a log-normal magnification PDF. The setup is otherwise the same as in 
Table HI 



•Scent 




Gaussian PDF 




0"68,GasB 
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1.74 
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1.09 


2 


0.279 
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2 


1.25 


1.33 


0.279 


6.47 X 10-3 


1.35 X 10-2 


1.04 


1.07 


3 


0.279 


9.50 X fO-^* 


1.94 X 10- 


2 


1.12 


1.11 


0.279 


8.78 X 10-3 


1.77 X 10-2 


1.03 


1.02 


4 


0.277 


1.22 X 10^2 


2.31 X 10" 


2 


1.07 


1.00 


0.274 


1.16 X 10-2 


2.34 X 10-2 


1.03 


1.01 



Table 5. Same as Tabled except that w is being constrained. 
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-1.003 


1.36 X 10-2 
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-2 
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1.57 


-0.991 
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1.75 X 10" 


-2 


1.01 


1.05 


2 


-1.006 


2.90 X 10-2 


6.03 X 10" 
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1.25 


1.27 


-1.000 


2.36 X 10-2 


4.94 X 10" 
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1.02 


1.05 
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-1.012 


4.89 X 10-2 


9.71 X 10- 


-2 
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1.10 
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-2 


1.02 


1.01 
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-1.000 
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-1 
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1.05 
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-1 
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1.00 



redshift bin is ~ 80% for a Gaussian PDF, but only ~ 10% 
for the log-normal PDF. This question has also been inves- 
tigated in the context of SNe la by DV06, who found that 
the variance is almost unchanged. We have confirmed their 
conclusion by running a similar Monte-Carlo simulation for 
a mock SNe la survey. The reason why the constraints on 
Sim from standard sirens, in contrast to standard candles, 
are obviously affected by the choice of the convergence PDF, 
is the large skewness (seen in Fig. [3| in the error distribu- 
tion of standard sirens (see more explicit demonstration of 
this point below). From Table [J] it is also clear that the 
improvement of the constraints on Sl,„ decreases with in- 
creasing source redshift. This feature could also be explained 
by the decrease in the skewness of the distribution with in- 
creasing source redshift. Here, the change is not driven by 
the intrinsic dispersion, but by the lensing dispersion. As 
we show in Figure [T] the convergence PDF is more skewed 
at low redshift, making the effect of non-Gaussianity more 
prominent. 

The conclusions for w are qualitatively very similar, ex- 
cept that the differences caused by the "wrong" PDFs are 
smaller, as shown in Table (5] In fact, the constraints with 
the log-normal PDF are almost identical to those from the 
"true" PDF, with a negligible bias, and with an overestimate 
of the 68% error by only 2%. We expect these qualitative 
conclusions, from the examples of Q,m and w, to remain valid 
for other parameters such as the Hubble constant h, or the 
dark energy equation of state evolution parameter dw/dz, 
whose impact is primarily through dhi^z), rather than the 
shape of the lensing PDF. 

Interestingly, as shown in Table |6] for ag, the erroneous 
use of the Gaussian PDF also degrades the constraints, while 
the log-normal PDF gives rise to better constraints, but 
causes a clear negative bias. At the lowest redshift, the bias 
is ~ 0.07, exceeding the forecast 68% constraint. The still 
larger, and positive, bias of for the Gaussian PDF seen in 
DV06 was not found in our calculations; instead, the peak is 
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Figure 7. The top panels show the contribution to § (a x2-likc 
quantity; see text for the definition) at z = 1 from events with 
different 5mtot (the deviation in the distance modulus from that 
in a homogeneous universe) . The middle and bottom panels show 
the change in the ^-contributions when the parameter is reduced 
or increased by one standard deviation creg.Gass- The left and 
right panels are for Cm and erg, respectively. 



of Q,m = 0.279 by at most « l/3rd of the 68% confidence 
level. This result is in agreement with DV06, and is not 
surprising, since the peak of the best-fit 0,^ distribution is 
determined primarily by the luminosity distance (or the ho- 
mogeneous distance modulus mo), rather than the shape of 
the 5m-distribution. The variance of the best-fit Jim, how- 
ever, becomes larger when the Gaussian or log-normal PDF 
is used in the likelihood analysis. The difference in the lowest 
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Table 6. Same as Table |4l except that erg is constrained and there are 400 events in each realization. 
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Figure 8. Distributions of the best-fit erg values, obtained by 
fitting the "true" magnification PDF with erroneous Gaussian 
PDFs, in two Monte-Carlo simulations that only differ in the size 
of the intrinsic scatter. The left (right) panel shows the results 
when the intrinsic scatter is typical of standard sirens (standard 
candles). Two thousand events are generated in each realization. 
Note that ag is negatively biased for standard sirens, while it is 
positively biased for SNe (see text and Figure[9]for an explanation 
of this result). 



slightly negatively biased in the lowest redshift bin (we find 
the mean, however, is almost unbiased). 

Investigations of the seeming inconsistency with DV06 
revealed, initially to our own surprise, that the sign of the 
(J8-bias is affected by the size of the intrinsic scatter, in 
addition to the skewness of the lensing PDF. To illustrate 
this point, in Figure |8] we plot the distributions of the best- 
fit (Tg from two Monte-Carlo simulations that only differ in 
the amount of the intrinsic scatter - the left panel assumes 
the intrinsic dispersion of a typical standard siren, while 
the right panel assumes that of a typical standard candle 
(cint = 0.1, as before). As the right panel of this figure 
shows, in addition to increasing the variance of the inferred 
erg distribution, a larger intrinsic dispersion also shifts the 
peak position to a higher value, in the case of the Gaussian 
PDF. This positive bias in the SN sample has a value of 
Acrg ^ 0.017. 

To understand why the erg-bias changes sign as the in- 
trinsic scatter is increased, we again consider the quantity ^ 
and its dependence on erg . In the bottom panels of Figure [9l 
we show (= ^(crg = 0.834) - ^(fiducial ag = 0.817), solid 
curves) for the case of the erroneous Gaussian PDF; the left 
and right panels are again for standard sirens and for SNe 
la, respectively. In agreement with the results of the Monte- 
Carlo simulations, AE from the integration of A^ is negative 
for SNe la, meaning the bias is positive, while AH for stan- 



dard sirens is positive, consistent with no bias. The top and 
bottom panels show Ptot{Smtot) and A{—2\nPlot{S'm[ot)), 
the two multiplicative factors comprising A^. For compari- 
son, we also show, by the dashed curves, the same quantities 
in a "true Gaussian" case, where the "true" convergence 
PDF (in addition to the PDF employed in the likelihood 
analysis) is assumed to be Gaussian, and therefore the esti- 
mate of (Tg should be unbiased. 

Compared with the Gaussian "true" distribution, the 
skewed "true" distribution produces a much higher proba- 
bility of strongly magnified events (solid curves in the top 
panels). These events contribute negative A^, i.e., they favor 
larger erg values (the solid curve is largely below the dashed 
curve at Snitot < 0). This is the effect that dominates in 
the case of SNe, as already argued by DV06. However, when 
the intrinsic dispersion is small, this negative contribution 
to H is, in fact, balanced by the decrease of the probability 
at 5mtot^0. In other words, the large erg value that is fa- 
vored by the events in the strongly magnified tail, is even 
more strongly disfavored by the lack of corresponding mildly 
demagnified events - leading to an overall negative bias in 
erg. In the SN case, when the intrinsic dispersion is large, 
there is no longer a lack of mildly demagnified events to be 
disfavored, because the expected number of these events is 
still dominated by the intrinsic scatter, and is insensitive to 
erg — the net effect, in this case, is the positive bias in erg, 
driven by the highly magnified events. 



4.5 Varying two cosmological parameters 

In our calculations above, we varied a single parameter, 
while all others were held fixed. To take into account the cor- 
relations among different cosmological parameters, we next 
consider simultaneous constraints on two cosmological pa- 
rameters. Ideally, all of the cosmological parameters should 
be allowed to vary simultaneously, for a full study of correla- 
tions. Unfortunately, with our simple Monte-Carlo method, 
this is computationally prohibitively expensive (and would 
require the use of Monte Carlo Markov Chains, or a similar 
technique). In this study, for simplicity, we instead chose to 
study only two typical correlations as illustrative examples. 
One of the correlations, between Qrn and w, is expected to 
be large, since as we have shown, both parameters draw their 
constraints from the mean distance modulus. The other cor- 
relation, between Jim and erg, is expected to be small, as the 
constraints on erg are purely from the cosmological depen- 
dence of the shape of the convergence PDF. We focus on the 
lowest redshift bin of 2cent = 1, where the assumed shape 
of the convergence PDF makes the largest difference. One 
hundred events are generated in each realization. 
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Table 7. The peak, 68% and 95% marginalized errors of and w for Gaussian and true PDFs. Each realization produces 100 events. 



Parameter 




Gaussian PDF 

''"eS.Gass 0"95,Gass 




True PDF 

""68, true 


fgS.truc 


0"68,tpuo 




w 


0.280 
-0.986 


1.44 X 10^2 
3.69 X 10-2 


2.87 X 10-2 
7.78 X 10-2 


0.277 
-1.000 


1.22 X 10-2 
2.92 X 10-2 


2.69 X 10-2 
6.18 X 10-2 


1.18 
1.26 


1.07 
1.26 



Table 8. The peak, 68% and 95% marginalized errors of Qm and trg for Gaussian and true PDFs. Each realization produces 100 events. 
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Figure 9. The top panels show the probability distributions of 
Sratot at ^ = 1. The solid and dashed curves assume that the 
true convergence PDF follows Pdo or a Gaussian, respectively. 
The middle panels show the change in —2\r\Pl^^{&m'^^^) when erg 
is increased from the fiducial value of 0.817 to 0.834 (the mean 
of erg from SNe la simulations when Gaussian PDF is used in the 
likelihood analysis). The bottom panels show the corresponding 
change in the x2— like quantity f . The left and right panels are for 
standard sirens and SNe la, respectively. 



The joint 2-D confidence contours and the marginal- 
ized 1-parameter distributions are plotted in Figure [10] (for 
simultaneous constraints on Q.m and w) and in Figure [11] 
(for simultaneous constraints on fim and erg). The values of 
the marginalized errors are also listed in Tables [7] and [8] 

Figure [TU] confirms that fi^ and w are strongly anti- 
correlated. As in the single parameter case, and w are 
both more tightly constrained with the true PDF, by 18% 
for fim and 26% for w (evaluated using their marginalized 
68% errors). The area enclosed by the 68% (95%) contour 
is larger in the Gaussian PDF case by 117% (92%). The 
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Figure 10. Upper panel: Contours enclosing 68% and 95% of the 
best-fit values on f!m and w in the lowest redshift bin for the 
true (solid line) and Gaussian (dashed line) PDF. In each case, 
10* realizations were produced of mock samples containing 100 
standard sirens, and Q,m and w were fit simultaneously, with all 
other parameters fixed at their fiducial values. Bottom left panel: 
Distribution of best— fit Om, marginalized over to. Bottom right 
panel: Distribution of best— fit ui, marginalized over Q,m- 



raw sensitivity to the combination of Q,m and w is therefore 
improved, by the presence of non-Gaussianity, by more than 
a factor of two. 

Figure [11] shows that the correlation between f2,„ and 
(jg depends on the shape of the convergence PDF. With 
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Table 9. The ratios of the 68% errors on the individual cosmo- 
logical parameters in the true and in the Gaussian cases, when 
the number of events in each realization is increased or decreased 
by a factor of 4 relative to our fiducial choices. 
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Figure 11. Same as Figure [Tol except that Om and erg are being 
simultaneously constrained. 

a Gaussian PDF, the parameters are almost uncorrelatecjfl, 
while with the true PDF, they are positively correlated. Here 
we offer a heuristic explanation for this positive correlation: 
with larger erg, the convergence dispersion is increased, mak- 
ing the convergence PDF more skewed, and its peak shift to 
lower k; to fit the same mtot, this requires a larger Jim to 
compensate. The marginalized error on Q.m is still smaller 
with the true PDF, but the difference of the 68% errors, 
only 19%, is considerably smaller than in the single parame- 
ter case. At the same time, the constraints on erg are slightly 
more degraded than those with a Gaussian PDF. Here, the 
area enclosed by the 68% (95%) contour is smaller in the 
Gaussian PDF case by 6% (10%). 



4.6 The importance of the number of events 

Another interesting question to ask is: how do our results 
- in particular, the tightening or degradation in the con- 
straints caused by non-Gaussianities - depend on the num- 
ber of standard siren events A'cvcnt? Depending on what 
statistic is being used, non-Gaussianities can become in- 
creasingly unimportant as A'^cvent — > oo. This is the case, for 
example, when one considers the distribution of the mean 
magnification {(5m) from Ncvant events (as is well known 
from the central limit theorem; see, e.g.. Figure 1 in Holz & 
Linder 2005). This statistic, however would be relevant only 
if one simply ignored lensing, and performed a least-squares 



^ In fact, Qm and erg are slightly anti-correlated, since both pa- 
rameters affect the lensing dispersion. 



Figure 12. A crude estimate for the number of SNe whose statis- 
tical constraining power is equivalent to that of a single standard 
siren. We have assumed that the magnitude errors of the two 
types of sources follow the same distribution, with the amplitude 
given by the convolution of intrinsic and lensing errors. 



fit to the observed Hubble diagram. This is different from 
the maximum-likelihood estimator we adopted here, which 
includes information beyond the mean magnification, and 
utilizes the full shape of the PDF. In fact, we expect that 
the MLE has a distribution for which the non-Gaussianity 
does not average away in the limit of a large number of 
events. 

In Table [9l we show the ratios of the 68% errors in the 
true and the Gaussian cases, when the number of events 
in each realization is increased and decreased, by a factor 
of 4, from the values adopted in § 14.11 As the table shows, 
the ratios are consistent with no dependence on %/ Acvcnt 
(while the absolute errors, which are not shown in the table, 
were found to scale approximately as V Acvcnt ) • Therefore, 
the general conclusion of this study, that the non-Gaussian 
shape non-trivially affects cosmological constraints, is valid 
independent of Acvont. This also demonstrates the impor- 
tance of choosing an optimal statistic: the information from 
the skewness in the PDF would be lost if only the mean 
distance modulus were used. 



4.7 The expected constraints from a LISA 
standard siren sample 

What is the actual expected error from LISA standard 
sirens on cosmological parameters, in the presence of lens- 
ing? Since we have not obtained these errors in the full 
multi-dimensional parameter space, we cannot fully answer 
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this question. However, the same question has been stud- 
ied exhaustively in the context of SNe - and, since the SNe 
constraints are from the same observable (the Hubble dia- 
gram), it makes sense to quote the expectations from stan- 
dard sirens relative to those from a SN sample. If the mag- 
nitude errors had the same distribution for these two types 
of sources, with standard deviations of CTss and (Tsn, then the 
statistical constraints from a single standard siren would be 
equivalent to that from A'an,! = {oaa / osu)'^ SNe. Assuming 
a constant (Jsn ~ 0.1 mag, and also that the lensing and in- 
trinsic errors add in quadrature, in Figure[l2]we show A^sn,!- 
As the intrinsic errors of standard sirens approach those of 
SNe, their relative advantage decreases monotonically with 
redshift, with iVen,i ^ 93 -5> 5.5 -)■ 2.5 at z 0.3 -s- 1 -)■ 1.7. 
(Peculiar velocities will limit the constraints from standard 
sirens below redshifts 2:;^0.3.) The main result of this pa- 
per could be stated by saying that a,t z — 1, these numbers 
for A'sn are increased by a factor of 1.5^ — 2.25 for flm, 
and by a factor of 1.8^ = 3.25 for w; whereas they are de- 
creEised by a factor of 0.8^ — 0.64 for ag. Note that while the 
actual numbers for A'sn.i are only crude estimates, these in- 
creases/decreases, which are due to the non-Gaussian shape 
of the lensing PDF, should be robust. At higher redshifts, 
these increases/decreases disappear. 

These estimates apply to a single redshift, and to a 
single cosmological parameter. The constraints from an ac- 
tual set of LISA standard sirens will, of course, depend on 
the redshift distribution of the sources. If most of the stan- 
dard sirens are at high redshifts (^^3; e.g., in Sesana 2004) 
then our results suggest that there will be no significant 
overall improvements from non-Gaussianity for the set of 
events If there are more events at low redshift, then non- 
Gaussianity will improve the constraints significantly. Hav- 
ing low-2 events remains a possibility (see, e.g., discussion of 
some of the uncertainties, and related references, in Lippai 
et al. 2009). HoUey-Bockelmann et al. (2010) recently found, 
using cosmological simulations of the BH merger trees, that 
most LISA-detectable merger events are below z = 2 (al- 
beit with high mass ratios, so their intrinsic di, error will be 
worse than for the equal-mass binaries adopted here). Fi- 
nally, even if most events are at ^^4, non-Gaussianities can 
still improve the overall constraints when the LISA data is 
combined with information from still higher redshift, such 
as cosmic microwave background (CMB) anisotropies. Since 
the CMB fixes the luminosity distance at 2 « 1,000, the 
lowest-redshift standard sirens will have the longest "lever- 
arm" and the contribution from these events to the com- 
bined LISA+CMB constraints will be enhanced significantly 
(a point emphasized in a different context by Hu & Haiman 
2003). 



5 CONCLUSIONS 

In this paper, we have shown that the sensitivity of low- 
redshift (z ~ 1) standard sirens to the parameters i},n and 
w are tightened, by a factor of 1.5 — 1.8, by the non-Gaussian 



' Of course, the fact that the standard siren sample extends to 
high redshift is still a useful complement to SN datasets, which 
are not expected to reach beyond z ~ 1.7 (Aldering et al. 2004). 



shape of the lensing magnification PDF, relative to a Gaus- 
sian PDF with the same variance. When these two parame- 
ters are constrained simultaneously, the improvement of the 
constraints, attributable to the skewness alone, is further 
enhanced, owing to the correlation between the parameters. 
We expect similar conclusions to hold for other parameters 
whose constraints are driven primarily by the changes they 
induce in the luminosity distance ^1,(2). Interestingly, the 
constraint on erg, which comes from the shape of the PDF 
itself, is, however, degraded by a factor of ~ 0.8. In our 
study, we relied on the shape of the lensing PDF described 
by the fitting formulae in Das & Ostriker (2006) . This shape 
is somewhat less skewed than the simulation results, suggest- 
ing that the improvement / degradation of the constraints 
has likely been underestimated. 

The improvements for flm and w are comparable to 
those that may be available from correlations with lensing 
measurements on larger scales, or from directly subtract- 
ing the lensing contribution of foreground galaxies. How- 
ever, these corrections from non-Gaussianities are "automat- 
ically" available, as long as the lensing PDF can be mod- 
eled ab-initio. We also found that at higher redshift, the 
effects from non-Gaussianity are less pronounced due to the 
reduced skewness of the lensing PDF. Overall, our results 
highlight the importance of obtaining an accurate and reli- 
able PDF of the point-source lensing magnification, in order 
to realize the full potential of standard sirens as cosmological 
probes. 
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